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1. Introduction 

Fuzzy spaces provide an attractive, alternative way to discretize spacetime and thus help 
regularize field theories on such spaces. Its interesting feature is the fact that this method 
preserves the symmetries of the continuum spacetimes even at the discrete level which are 
broken in lattice regularization techniques. 

Such fuzzy spaces are studied using the noncommuting algebra of functions defined 
on them. One such widely studied model is the fuzzy sphere Sp ^. Several works on 
numerical simulations of scalar fields and gauge fields on Sp have been done. 0, ^, ^, |7[ § |^ . 

There have also been interesting attempts to extend physics to other such fuzzy spaces 



with higher genuses [10] and more exotic looking surfaces |11, 12 1. But here we will consider 
only Sp. 

In Connes' approach to noncommutative geometry [^], the Dirac operator gains funda- 
mental significance as part of the spectral triple in formulating the spectral action principle. 
Such a spectral action has been considered recently |14| for the Dirac operator of a spin 
^ particle. The operator used by them corresponds to the one constructed in [^]. In this 
work we will consider Dirac operators constructed with the help of Ginsparg-Wilson(GW) 
algebras. This approach provides an elegant way to extend the construction of Dirac oper- 
ators to all spins as studied in Q, where it was also found that several such Dirac operators 
exist for the case of each spin j. 



We consider the spin 1 case and numerically compute the spectrum of the 3 Dirac 
operators. We also analytically compute the traces of these 3 operators and find that two 
of them have non-zero trace showing the existence of unpaired eigenstates(zero modes). 
One of them is traceless. Rather surprisingly, the spectrums are found to be different, both 
numerically and, analytically from the trace formulas, thereby establishing the unitary 
inequivalence of the three operators. The spectrum of the traceless Dirac operator is 
studied in detail. The eigenvalues are plotted as a function of their degeneracy. This is 
done for values of the cut-off, L for which we could compute the eigenvalues numerically. 
We then fit this data with curves for each L and find that a quadratic function fits it. This 
fit extends rather amazingly to all other values of L once we change the parameters in the 
fit which are functions of L. Thus this gives an universal relation between the energy and 
their degeneracy and helps predict the eigenvalues for any large value of L. 

The paper is organized as follows. Section 2 briefly describes the noncommutative 
algebra on S'^. We then recall the construction of Dirac operators using GW algebras 
in section 3. Here we also write down the 3 spin 1 Dirac operators which we will work 
with. The trace of these operators are found analytically. Their inequivalence provides a 
simple proof for the unitary inequivalence of the three operators. In section 4 we recall 
the spectrum of the spin ^ Dirac operator and we study the properties of the spectrum of 
the spin 1 Dirac operator. In particular we provide an elegant way to count the number 
of zero modes of the traceless spin 1 Dirac operator for each cut-off L. The numerical 
results are presented in section 5. We conclude in section 6 with a few remarks and further 
speculations. 



2. Geometry of S% 

The algebra for the fuzzy sphere is characterized by a cut-off angular momentum L and is 
the fuU matrix algebra Mat{2L + 1) = M2L+1 of (2L + 1) x (2L + 1) matrices. They can 
be generated by the (2L -|- l)-dimensional irreducible representation (IRR) of SU{2) with 
the standard angular momentum basis. The latter is represented by the angular momenta 
Lf acting on the left on Mat{2L + 1): If a G Mat{2L + 1), 

Lfa = Lia (2.1) 

[LtLf]=i^ijkLj: (2.2) 

(Lf)2 = L(L + l)l (2.3) 

where Lj are the standard angular momentum matrices for angular momentum L. 
We can also define right angular momenta Lf^: 

Lfa = aLi,ae M2L+1 (2.4) 

[L^,Lf] = -ieijkL^ (2.5) 
(Lf)2 = L(L + 1)1 (2.6) 



We also have 

[Lf,Lf]=0. (2.7) 

The operator Li = Lf — is the fuzzy version of orbital angular momentum. They 
satisfy the SU{2) angular momentum algebra 

[Ci,Cj] = ieijkCk (2.8) 

In the continuum, 5^ can be described by the unit vector x G S'^, where x.x = 1. Its 
analogue on Sp is or such that 

lim ^— = X,. (2.9) 

This shows that if do not have continuum limits. But £j = Lf — does and becomes 
the orbital angular momentum as L — > oo: 

lim Lf - Lf = A (2.10) 

L^oo 

3. Construction of the Dirac Operators 

In algebraic terms, the GW algebra A is the unital * algebra over C , generated by two 
^-invariant involutions r,r'. 

^ = {r, r' : = r'2 = 1 , r* = r , r'* = r'} (3.1) 

In any * -representation on a Hilbert space, * becomes the adjoint f. 
Consider the following two elements constructed out of r,r': 

ri = ^(r + r'), (3.2) 
r2 = ^(r-rO. (3.3) 

It follows from Eq.(^) that {ri,r2} = 0. This suggests that for suitable choices of F, 
r', one of these operators may serve as the Dirac operator and the other as the chirality 
operator provided they have the right continuum limits after suitable scaling. 

For the spin 1 case the combination which leads to the desired Dirac and chirality 
operators were found in Q and they are 

D, = L , (3.4) 

D2 = L I^Ei^l-I^^ (3.5) 

and 

L rR 



Dz = L[ ^ . (3.6) 



with 



71 = ^ , (3.7) 



72 = ^ j (3.8) 

and 

73 = {^^^) (3.9) 
as their corresponding chiraHty operators. In the above equations 

r 2(S.L^ + L + 1)(S.L^ + 1) - (L + 1)(2L + 1) 



(3.10) 



(L + 1)(2L + 1) 

^ _ 2(-g.L^ + L + l)(-g.L« + 1) - (L + 1)(2L + 1) 
^ " (L + 1)(2L + 1) ' ^"^-^^^ 



,i _ 2(S.L^-L)(S.L^ + 1)-L(2L + 1) 



= L(2L + 1) (^-1^) 



^ ^ 2(g.L^ + L)(g.L^ - 1) - L{2L + 1) 
L(2L + 1) 



(3.13) 



r _ -2(S.L^ - L)(S.L^ + L + 1) - L(L + 1) 



and 



(3.15) 



^ ^ 2(S.L^ + L)(-S.L^ + L + 1) - L(L + 1) 
^" L(L + 1) 

The operators in Eg. ( ^. 10 )-Eq. dsTlSl ) are generators of GW algebras and are obtained from 
left and right projectors to eigenspaces of the total angular momentum, L + S, where S 
are the matrices representing the spin 1 representation of SU{2). 
The continuum limits of Eq.(|3.4|)-Eq.(|3.6|) are 



Di = - (S.x)^ + 2) + 2(S.x) + S.x}, (3.16) 

D2 = {t.C - {t.xf + 2) - 2(S.x) - S.x} (3.17) 

and 

= (S.x)^ + 2. (3.18) 
The corresponding chirality operators in the continuum are 

71 = (S.x)2 + (S.x)-l, (3.19) 

72 = (S.x)2 - (S.x) - 1 (3.20) 

and 

73 = l-2(S.x)2 (3 21) 

respectively. 



Dirac Operator 




Let 


Di 


4L(2L + 1) 


2L(5L + 1) 


D2 


-4L(2L + 1) 


-2L(5L + 1) 


D3 









Table 1: Traces of the 3 Dirac Operators 



The trace of the Dirac operators 

The trace of the Dirac operators in Eq.(|3.4D-Eq.(|3.6D can be computed analytically by using 
the formula 

tr{A(^ B) = tr{A).tr{B) (3.22) 

where A and B are square matrices. Since the Dirac operators we construct act on 
Mat{2L + 1) (8) C^, they are of the form of tensor products and hence we can apply this 
formula to analytically compute their traces. 

The trace is a rotationally invariant object leading to 

tr{{Lff)=tr{{L^f)=tr{{L^f) (3.23) 

and 

tr{T.l) = tr{T,l) = tr(S|) = 2. (3.24) 

The above equations hold due to the fact that the three generators of any representation 
of the SU{2) algebra have the same trace because of rotational invariance. 

The trace of [L^)"^ varies according to whether L is integer or half- integer. When L is 
an integer 

tr{{Lff) = \L{L + l){2L + lf (3.25) 
and when L is an half-integer 

tr{{Llf) = hiL + 1)(L + 2)(2L + 1). (3.26) 

The same formulas hold when the left operators in the above equations are replaced by 
right operators. It is simple to see that Yl.L^ and Yl.L^ are traceless. Using these identities 
we write down the traces of our 3 Dirac operators in Table(2) 

The trace of the Dirac operator is the sum of its eigenvalues. The availability of 
these exact trace formulas are helpful in verifying the spectrum of these operators found 
numerically. 

The operators Di and D2 have non-zero trace implying the existence of unpaired 
eigenstates or zero modes. 

To check the unitary equivalence of the 3 Dirac operators, it is a necessary, though 
not sufficient condition that the traces of the 3 operators be the same. Since the trace 
formulas show the traces are not the same, they provide an analytic proof for the unitary 
inequivalence of the 3 Dirac operators confirming numerical results. 



4. Analytic results of the spectrums of the spin | and spin 1 Dirac oper- 
ators 



The spectrum of the spin | Dirac operator can be found analytically [^]. In the GW 
approach to constructing the Dirac operator, the spin i system has the same spectrum both 
in the continuum and the fuzzy level. To illustrate the method of finding the spectrum, we 
consider the spin ^ Dirac operator in the continuum: 

Di=a.C + l. (4.1) 

2 

In the above equation C is the orbital angular momentum got by taking the continuum 
limit of — L^. a are the spin ^ Pauli matrices. The total angular momentum J given 
by 

commutes with the Dirac operator. We can use its eigenvalues to label the eigenstates of 
the Dirac operator. For given cut-off angular momentum L, the spectrum of the orbital 
angular momentum is given by 

£g {0,1,--- ,2L}. (4.2) 
Given this we can find the spectrum of the total angular momentum J to be 

,2L-i,2L + i}. (4.3) 

Each value of the total angular momentum J can be got from two different orbital angular 
momentum except the top mode whose J value is 2L + |. From this we can count the total 
number of eigenvalues for a given cut-off L with the help of the following sum: 

2(2i)+2L + - = 2(2L + l)2. (4.4) 



The spectrum of the Dirac operator in Eq.(4.1) can be got by noting that this operator 
can be written as 

Di=J^-C^ + \. (4.5) 

2 4 

As J'^,C'^ = 0, we can write the spectrum of Di as 

J 2 

1 

Spectrum of Di = j{j + 1) - l{l + 1) + - . (4.6) 

2 4 

As mentioned before each j comes from two different I values except the top mode. Thus 
we have for the spectrum of Di: 



The spectrum has the chiral nature as expected. Note that there are no zero modes for 
the spin ^ Dirac operator. The computation of the spectrum in the spin ^ is easy due 
to the form of Di as given by Eq.(|4.1|). This however is not true for the Dirac operator 

of the spin 1 case given by Eq.( |3.18| ). This is due to the presence of the term S.x which 
does not commute with Ti.C making the analytic computation difficult. This is the reason 
why we take to numerical methods to achieve this. Nevertheless we can still get some vital 
information about the spectrum of the spin 1 Dirac operator by analytic methods. 
The total angular momentum J given by 

J = S + £ 



commutes with the Dirac operator in Eq.( |3.18 ) just as the corresponding total angular 
momentum does in the spin ^ case. The spectrum of the orbital angular momentum C is 
the same as in the spin ^ case given by Eq.( [4.2D . The spectrum of J is now given by 

J e {0, 1, 2, • • • , 2L - 1, 2L, 2L + 1}. (4.8) 

In this case each value of J comes from three different orbital angular momenta namely 
j — 1, j and j + 1 except three j values, j = comes from only one state, j = 2L comes 
from 2 states and j = 2L + 1 comes from only one state. These are easy to check as they 
involve the simple angular momentum addition rules. With this information we can count 
the number of eigenvalues for each cut-off L with the following sum: 

j=2L-l 

1+ 3(2j + l) + 2(4L + l) + 2(2L + l) + l = 3(2L + l)2. (4.9) 

i=i 

This is exactly the number of eigenvalues we expect from each cut-off L for the spin 1 case 
as this is the size of the matrix for the Dirac operator for each L. These arguments can be 
easily extended to the Dirac operators of all spins but we will not do so here. 

4.1 Number of positive eigenvalues and Zero modes for the spin 1 Dirac oper- 
ator 

Out of the three Dirac operators in the spin 1 case we will consider the traceless Dirac 
operator (See Table ||). The trace equation gives us an easy and elegant way to count the 
number of different non-zero positive and negative eigenvalues as well as the number of 
zero modes for each cut-off angular momentum L. 

The zero modes can be counted as follows: j = comes from just one orbital angular 
momentum state and so it cant result in a positive or negative eigenvalue of and hence 
it must only be due to the traceless nature of the Dirac operator. This contributes 1 zero 
mode for each L. Similar argument holds for j = 2L -|- 1 which contributes 2(2L -|- 1) -|- 1 
zero modes for each L. For values of j between 1 and 2L — 1 there is a contribution of 
2j + 1 zero modes for each of the j values. Summing all this we find that there are exactly 
(2L -|- 1)^ -|- 2 zero modes for each L. 

In a similar way we can find the number of positive eigenvalues. When we do this we 
find there is a contribution of 2j-|-l eigenvalues for values of j between 1 and 2L. Summing 



these we get 4L^ + 4L. As the Dirac operator is traceless, the same argument holds for 
the negative eigenvalues giving a total of + 4L eigenvalues for each L. It is easy to see 
that the sum of the positive, negative eigenvalues and zero eigenvalues give 3(2L + 1)^ as 
the total number of eigenvalues as expected for each cut-off L. 

These arguments can again be easily extended to the spectrum of higher spin Dirac 
operators on Sp. It is also easy to see that there are no zero modes for half- integral spin 
systems on Sp as none of the Dirac operators for half-integral spin systems are traceless. 
We will not discuss them any further in this work except for a few remarks in the end. 

Finally we count the number of different positive eigenvalues we expect to find for the 
spin 1 Dirac operator for each cut-off L. Since there are 2L + 2 values the total angular 
momentum j can take, out of which 2 of them can only contribute to the zero modes for 
each L, we can conclude that there are 2L different positive eigenvalues for each L. The 
degeneracies of each of them can easily be read off as 2j + 1 according to the corresponding 
value j takes. 

4.2 Remarks on the other two Dirac operators for the spin 1 case 



So far the arguments in this section were for the traceless Dirac operator in Eq.( 3.1q ). 
These arguments do not hold for the Dirac operators in Eq. (|3.16| ) and Eq. (|3.17j ) as they 
have positive and negative traces respectively. These are given in table ||. 

Consider the Dirac operator with the positive trace whose continuum value is given by 
Eq.( p.l6| ). In this case too we have for the spectrum of the total angular momentum 

Spec Jg {0,1,- •• ,2L-l,2L,2L + l} (4.10) 

as before. However in this case we cannot say that the states corresponding to j = 2L + 1 
and J = correspond to zero modes. This is because of the non-zero trace. They now have 
some positive energy say Eq and E2l+i- We then have the following equation 

+ (4L + 3)S2L+i = 4L(2L + 1) (4.11) 

for integral values of L and 

^0 + (4L + 3)^2L+i = 2L(5L + 1) (4.12) 

for half-integral values of L. The 4L -|- 3 states with energy E2L+1 correspond to unpaired 
eigenstates. If |£'2l+i) is the state with energy E2L+1, then these states will be of the form 
7^'^|£'2L-f 1) where 7 is the chirality operator given in Eq.( [3.19 ) and k is an integer. In the 



spin ^ case these states, with the corresponding chirality operator for the spin ^ system, 
will equal itself as 7^ = 1 for the spin ^ case. It can be easily seen from Eq.( |3.19| ) that 
this is not true for the spin 1 case. So we get the possibility for a number of states with the 
same energy. With this note, we analyze only the spectrum of the traceless Dirac operator 
in what follows. 

Having studied the general nature of the spectrum for the spin 1 Dirac operator, we 
compute the eigenvalues numerically in the next section. In particular we will find a 
relation between the eigenvalue and its multiplicity for a given cut-off L. This is equivalent 
to finding the eigenvalues as a function of total angular momentum j for each cut-off L. 
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Figure 1: Plot of the energy eigenvalues along with the best fit curves for iV = 21, = 35 and 
N = 45. 

5. Numerical Results 

We compute the eigenvalues of the the three Dirac operators in Eq.( |3.4D -Eq.( |3.6D numeri- 
cally. The size of each of these operators is 3A^^ where N = 2L + 1 and L is the cutoff. It is 
clear from the dimensions of these matrices (~ 9A^^) that we cannot go to arbitrarily large 
values of A^. Even for L = 22, the size of the matrix becomes 6075 x 6075 which is difficult 
to handle numerically within the resources available to us. For large values of N number of 
computational steps increase which will lead to growth of systematic error. However, the 
patterns emerging from the spectrum we computed so far strongly suggest what the behav- 
ior would be at higher values of N. This circumvents computational problems and helps 
us predict the behavior as we go close to the continuum. This is particularly important 
given the problems in handling very large matrices. 

The nature of the spectrum was discussed in the previous section and we confirm those 
results numerically. The spectrum of is similar to that of fermions with equal number 
of positive and negative energy eigenvalues. This is a reflection of the existence of the 



chirality operator given by Eq.(3.9), which anticommutes with D3. Apart from the non- 
zero eigenvalues there also exist a number of zero modes. We find exactly (2L + 1)^ + 2 zero 
modes for each cut-off L as we explained in the previous section. The number of positive 
eigenvalues is also as expected. 

We work with only the positive eigenvalues of the spectrum. As the operator is traceless 
we have the same pattern for the negative eigenvalues and so we do not use them to fit 
curves. Then we find the degeneracies of each of the positive eigenvalues. Note from the 
discussion in the earlier section, that there can only be odd degeneracies for our system 
as the total angular momentum j takes integral values. The plot for the energy vs the 
degeneracies is shown in figure 1. It shows the data points for three different values of N 
(namely = 21, = 35 and = 45) along with the best fit curves. 

By inspection we found the curve has a mirror symmetry about some principal axis. 



Next we try to find a universal curve that will fit the data(eigenvalues) of different 
cutoffs, just by changing the value of the cutoff. To this end we analyzed the data in a 
rotated frame in which the data was found to have reflection symmetry around the rotated 
y-axis. Given that the data for small {E,g) is independent of the cutoff(this can be seen 
in figure 1 where for small values of g the three sets of data points lie almost on top of 
each other on a straight line), we found a unique rotation angle to rotate all the results for 
different cut-off L. After observing the reflection symmetry of {E',g') we tried to fit the 
data with a polynomial with only even powers. To our surprise we found an excellent fit 
with just a parabola for all different cutoff values. Higher powers in the function did not 
make any further improvement in the fitting. The parameters of the parabola run with the 
cut off. We also find excellent fit for these parameters as a function of the cutoff L. 

We now elaborate this method. The plot of {E,g) is rotated to a new set of variables 
{E',g'). This set of points is then fitted with the curve 

E' = a{g' + rjf + f3. (5.1) 

Here a, /3 and rj are expected to vary with the cut-off L. The relation between (E' , g') and 
{E, g) is given by 

' E'\ ( cose sme\ ( E , . ^. 

(5.2) 



g' I \ — sin 9 cos t> I \ g 




The angle 9 = 2.26159radians. This angle is a constant for different values of L. This can 
be seen as a consequence of the data points lying on top of each other for small values of 
g as seen in figure 1. Note that E' and g' are not energies and degeneracies respectively. 



We just need to use the transformation in Eq.(5.2) to get the relation between the energies 
and the degeneracies. 

Our next task is to find a, /3 and tj as functions of = 2L + 1. They are found by 
fitting the quadratic form (Eq.(^?l|)) to the rotated curves for different values of A^. We 
find them to follow simple relations. These are shown in the figures 2. The exact functions 
we found were: 

0.863569 , , 

« = WMm^ - 0.00141635, (5.3) 

/3 = -1.45123A^ + 0.333497 (5.4) 

and 

r] = -1.16288Ar - 0.555529. (5.5) 

The numbers may look uninteresting but if we could fit these functions after we find these 
parameters for more values of A^ we could converge onto some special numbers. We did 
not attempt this in this work. The relations are simple enough to imply something deeper 
in the spectrum. More exact numbers could help in the quest for an analytic solution of 
this problem. 

We can now write down the exact relation between E and g based on our numerical 
fits: 
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Figure 2: Parameters a, j3 and 77 as a function of N . 



where 

a = acos^ 6, (5-7) 
h{g) = 2acos6 {smOg + T]) + sm.9 (5.8) 

and 

c{g) = a (sin 9g + 77)^ + /5 — cos 9g. (5-9) 

Having found this relation between E and g, we can now find the eigenvalues for 
arbitrarily large values of N . If we diagonalize -D3 for such large values of N it would take 
a lot of memory on the computer and is subject to a lot of numerical error. But we can 
get around this with our relation between E and g. Figure 3 shows the eigenvalues for 
N = 60. Note that though the curve looks continuous, we have seen in the previous section 
that degeneracies are allowed to take only odd integral values. The maximum degeneracy 
for a given is 2N — 1. Starting from 3 we can allow g to vary till 2N — 1 through odd 
integers and find the corresponding eigenvalues using Eq.( |5.6| ). In an equivalent manner we 
can find the energy eigenvalues as a function of the total angular momentum j by simply 
substituting g = 2j + 1 m Eq. (|5.6|) . 

6. Conclusions 

The spectrums of the spin 1 Dirac operators are found numerically. The three operators do 
not have the same spectrum making them unitarily inequivalent. This may have interesting 
consequences which we plan to explore in the future. The fermionic character of the 
spectrum is noteworthy as there exists no such higher dimensional analog in the Minkowski 
case. We expect this behavior also for higher spin Dirac operators on S'p as they all come 
paired with an anticommuting chirality operator. 

The universal relation between the energy eigenvalues and their degeneracies we find 
in Eq. (|5.6| ) is indeed remarkable as it circumvents what would have been an almost im- 
possible computational problem involving large matrices. This now allows us to find the 
eigenvalues for any arbitrarily large cut-off L. The simple relation between the eigenvalues 
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Figure 3: Plot predicting the energy eigenvalues for = 60. 

and their degeneracies seem to suggest some connection with the underlying symmetry in 
the problem. 

Having obtained the spectrum for the spin 1 Dirac operator for arbitrarily large cut- 
off L, we can go on to find the partition function for a system of particles occupying 
these energy levels. Assuming fermionic statistics we can compute several thermodynamic 
quantities for this system. Several interesting features arise and these are reported in J15| . 

A quantum particle on the continuum sphere, S"^ has energy eigenvalues given by 
/(/ + 1). These are the eigenvalues of the Laplacian on the sphere which is a second 
order differential operator. The eigenvalues of the square of the continuum limit of the 
spin i Dirac operator on Sp [|l6|, ||, ^, 17] also gives a spectrum similar to that of the 
standard Laplacian on S'^ apart from a additional constant. This additional constant can 
be interpreted as the scalar curvature according to the Lichnerowicz formula for the square 
of a general Dirac operator. In the Minkowskian case this is analogous to the square of 
the Dirac operator giving the Laplacian on that space. This leads to each component of 
the Dirac spinor satisfying the Klein-Gordon equation. We can view the Laplacian of the 
standard sphere as an analog of the Klein-Gordon equation on the sphere as this gives the 
SU{2) covariant dispersion relation on S'^. Note that we can add additional constants to 
this Laplacian as they are rotationally invariant. This however is not true for any of the 
3 spin 1 Dirac operators on Sp as their continuum limits, given by Eq.( p.l6D -Eq. (|3.18|) , 
contain S.x terms which makes the square of these operators look complicated. (Note 
that we do still get l{l + 1), but with additional terms containing S.x which makes the 
analytical computation of the spectrum difficult.) Thus the spectrum of their squares are 
not the standard one making the study of these deviations interesting as there exist no 
counterparts on higher dimensional Minkowskian space. 

We are also computing the spectral action of these Dirac operators. This will be 
compared with the spectral action of the spin ^ Dirac operator on the continuum sphere. 
This kind of analysis was carried out recently where interesting connections were made 



o 
o 

o 
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with cosmology. The results will be reported in a future work 

It was found in Q, that for a given chirality operator there exist several different Dirac 
operators. However this was done in the continuum limit and we have not found their fuzzy 
analogs. These will most certainly not be unitarily equivalent. This seems to be a new 
property of spin systems on Sp and they need to be studied further. 
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